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I.  INTRODUCTION 


The  chemical  reaction 

0  +  H2(v,j)  -  OH(v',j')  +  H  (1) 

and  the  subsequent  radiative  decay  of  the  v'  »  1  vibrational  mode  of  the  OH 
molecule  at  2.8  pm,  contribute  to  the  infrared  signature  of  certain  high- 
altitude  missile  plumes.  This  reaction  occurs  when  ^  molecules  present  in 
the  plume  collide  at  high  velocity  with  ambient  atmospheric  0  atoms.  Model 
calculations  of  the  radiation  generated  by  this  process  require  as  an  essen¬ 
tial  input  the  state-to-state  cross  sections,  Or(v,  j;  v'  -  l,j"),  of  this 
reaction.  The  range  of  the  initial  vibration  and  rotation  quantum  numbers 
v  and  j,  and  the  relative  translational  energy  E  of  the  reactants  for  which 
these  cross  sections  are  required,  are  determined  by  the  velocity  and  tempera¬ 
ture  of  the  plume. 

The  velocity  of  the  plume  is  the  vector  sum  of  the  vehicle  velocity  and 
the  rocket  exhaust  velocity.  This  will  be  a  maximum  in  the  retro-firing 
orientation  where  the  two  component  velocities  add.  At  high  altitudes  we  can 
assume  a  vehicle  velocity  of  about  ~6  km/sec  and  a  rocket  exhaust  velocity  of 
about  ~3  km/sec.  In  addition,  the  random  thermal  motion  of  the  reactants  will 
contribute  ~1  km/sec,  so  the  maximum  relative  collision  velocity  between  the 
H2  molecules  in  the  plume  and  the  atmospheric  0  atoms  will  be  of  the  order  of 
~10  km/sec.  This  corresponds  to  ~22  keal/mol  of  relative  translational  energy 
between  the  (0,^)  reactants. 

The  temperature  of  the  plume  decreases  rapidly  due  to  expansion  after 
leaving  the  rocket  nozzle.  Along  a  given  streamline  the  collision  frequency 
decreases  and  at  some  point  is  not  sufficient  to  maintain  thermal  equilibrium 
between  the  vibrational,  rotational  and  translational  modes.  The  vibrational 
modes  are  the  first  to  decouple  and  are  said  to  freeze  out  at  a  vibrational 
temperature  Ty.  Further  along  the  streamline  and  at  a  lower  temperature  TR, 
the  rotational  modes  freeze  out.  Typically,  the  rotational , temperature  will 
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be  quite  low.  We  have  assumed  that  T  <  200°K.  At  this  temperature  only  the 
j  *0,1  rotational  states  of  H2  will  be  significantly  populated.  The 
vibrational  temperature  Ty  is  higher  than  the  rotational  temperature,  but  is 
certainly  much  less  than  the  temperature  inside  the  rocket  engine,  which  is 
about  ~3000®K.  Even  using  this  very  crude  upper  limit,  the  conclusion  is  that 
only  the  v  *  0  and  v  *  1  vibrational  modes  of  H2  will  be  sufficiently 

populated  to  contribute  significantly  to  the  reactive  cross  section. 

We  report  here  the  results  of  a  Monte  Carlo  classical  trajectory  study  of 
the  state-to-state  cross  sections  for  reaction  (1).  These  cross  sections  are 
determined  in  the  range  of  relative  kinetic  energy,  from  the  lowest  threshold 
of  4  kcal/mol  to  an  upper  value  of  25  kcal/mol  and  for  initial  vibration- 
rotation  quantum  numbers  in  the  range  v  ■  0,1;  j  ■  0,1.  The  potential  surface 
used  in  these  calculations  is  discussed  in  Section  II.  The  calculational 
procedure  and  the  cross  section  results  are  given  in  Section  III. 


II.  POTENTIAL  ENERGY  SURFACE 


The  potential  function  used  in  this  study  is  the  same  potential  employed 
by  Johnson  and  Winter  (JW)^  in  a  previous  classical  trajectory  study  of  the 
(0,^)  system.  This  function  is  a  London-Eyring-Polanyi-Sato  (LEPS)  surface 
with  a  single  Sato  parameter  adjusted  by  trial  and  error  to  yield  a  computed 
rate  constant  in  agreement  with  experimental  data  around  320  K.  Rate  con¬ 
stants  computed  with  this  potential  were  in  good  agreement  with  the  best 
experimental  measurements  over  the  temperature  range  300  K  to  1000  K.  This 
potential  was  also  used  to  calculate  the  rate  constant  for  the  reaction  0  + 

»  1)  +  OH  +  H,  where  the  H2  ^as  a  single  quantum  of  vibrational 
l  2 

energy.  Light  has  since  measured  this  rate  at  T  *  302  K.  The  computed 
result  is  in  reasonably  good  agreement  with  the  experimental  value. 

O 

Recently,  two  groups,  Howard,  McLean  and  Lester  (HML)  ,  and  Walsh, 
Banning,  Raffenetti  and  Bobrowicz  (WDRB)^,  have  published  ab  initio  configur¬ 
ation  interaction  calculations  of  the  potential  surface  for  the  0(^P)  + 
H2(1Eg+)  reaction  which  can  be  compared  to  the  JW  LEPS  surface.  The  saddle 
point  properties  of  the  three  surfaces  are  compared  in  Table  1.  We  see  that 
the  JW.  potential  and  the  WDRB  potential  both  have  the  same  barrier  height  of 
12.5  kcal/mol,  while  the  HML  potential  is  14.0  kcal/mol.  Shinke  and  Lester 
have  fit  the  HML  surface  to  a  56-parameter  function  by  least  squares  and  used 
it  in  a  classical  trajectory  study'’  of  the  0  +  H2  system.  Their  calculated 
rate  constants  k(v)  for  the  reaction  0  +  H2(v)  +•  OH  +  H  are  less  than  the 
experimental  values  over  the  entire  temperature  range  from  300  K  to  1000  K  for 
the  case  v  -  0  and  are  also  less  than  Light's  experimental  value  for  the  case 
v  ■  1.  This  strongly  indicates  that  the  14.0  kcal/mol  barrier  height  of  the 
HML  potential  is  too  high. 

For  noncollinear  geometries  the  potential  surface  splits  into  two  dis- 
3  3 

tinct  surfaces  of  A"  and  k"  symmetry.  WDRB  have  calculated  both  these 

3 

surfaces.  In  the  saddle  point  region,  the  only  different  property  of  the  k' 
and  k"  surfaces  is  the  bending  potential  curve.  The  bending  curve  of  the 
JW  surface  is  approximately  the  average  of  these  two  curves^  and  thus 
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Table  1.  Comparison  of  potential  surfaces  for 
0  +  H2  in  the  saddle  point  region. 


HE 

rHH 

r0H 

kcal/mol 

A 

A 

Johnson,  Winter 
LEPS 

12.5 

0.95 

1.12 

Howard,  McLean 
and  Lester, 
ab  initio 

13.95 

0.95 

1.15 

Walch,  Dunning, 

12.5 

0.93 

1.21 

Raffenetti  and 
Bobrowicz , 
ab  Initio 


represents  a  very  good  single-surface  compromise  to  using  two  surfaces  for 
dynamical  calculations. 

Given  these  comparisons  of  the  JW  surface  with  the  ab  initio  surfaces  and 
the  good  agreement  of  the  JW  calculated  rate  constants  with  the  experimental 
values,  and  taking  account  of  the  uncertainties  in  both  the  ab  initio  poten¬ 
tials  and  the'  measured  rate  constants,  we  believe  that  the  simple  JW  LEPS 
potential  is  a  least  as  good  as  the  ab  initio  potentials  and  is  a  good 
approximate  surface  for  use  in  the  present  classical  trajectory  calculations. 


III.  CALCULATIONS  AND  RESULTS 


A.  TOTAL  REACTIVE  CROSS  SECTIONS 

Muckermann's  classical  trajectory  program,  CLASTR,^  which  is  based  on  the 
analysis  of  Karplus,  Porter  and  Sharma^,  was  used  to  calculate  the  cross  sec¬ 
tions.  Essential  inputs  to  the  program,  necessary  to  define  the  initial  state 
of  the  reactants,  are  the  separation  of  the  reactants  rQ,  the  relative  kinetic 
energy  E,  and  the  vibration-rotation  quantum  numbers  v  and  j.  Two  other  re¬ 
quired  input  quantities  are  bma y,  the  maximum  value  of  the  impact  parameter, 
and  N-j-qt*  t*ie  total  number  of  trajectories  to  be  calculated.  The  program 
computes  NT0T  trajectories,  selecting  initial  conditions  by  the  quasi- 
classical  Monte  Carlo  procedure;^  it  classifies  each  trajectory  as  reactive  or 
not  reactive;  and  it  calculates,  by  the  histogram  binning  method,  the 
vibration-rotation  quantum  numbers  v"  and  j"  of  the  product  diatomic 
molecule.  The  output  of  the  program  is  N^(E;  v,j;  v*,j*)f  which  is  the  number 
of  reactive  trajectories  out  of  the  total  number  NTQT  (E,v,j)  that  started  in 
state  v,j  and  ended  in  the  product  state  v",j"’. 

A  series  of  calculations  on  the  (O,^)  system  were  carried  out  using  this 
program.  The  initial  separation  in  all  cases  was  rQ  *  4  A,  the  kinetic  energy 
was  sampled  in  the  range  4  <  E  <  25  kcal/mol,  and  the  initial  vibration- 
rotation  quantum  numbers  were  v  ■  0,1  and  j  *  0,1. 

For  maximum  computational  efficiency,  bmax  should  not  be  any  greater  than 
the  maximum  impact  parameter  for  which  a  reaction  occurs.  This  varies,  de¬ 
pending  on  E,  v  and  j.  We  were  able  to  make  use  of  our  previous  calculations 
on  this  system*  to  obtain  reliable  estimates  of  bMX.  For  v  ■  0,  b^x  was 
varied  between  0.6  A  and  1.4  A  and  for  v»  l  it  was  varied  between  1.2  A  and 
1.9  A  depending  on  the  energy  E.  was  also  varied,  depending  on  the 

values  of  E  and  v.  For  v  -  1,  N^,q^,  ■  (400  -  1200)  trajectories,  and  for 
v  »  0,  Ntqt  *  (1200  -  2000)  trajectories.  In  general  the  larger  number  of 
trajectories  were  used  at  low  energies,  where  the  reactive  cross  section  is 
small,  in  order  to  maintain  enough  reactive  trajectories  to  calculate  with 
reasonable  statistical  accuracy  the  final  vibration-rotation  distribution. 
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The  toCal  reactive  cross  sections  for  a  specific  collision  energy  E  and 
initial  vibration-rotation  state  v,j  is  computed  by  the  formula^ 


°r(E.v.  j) 


,b„<  VE>v'«/!WE’"-i) 


(2) 


where 


Nr(E,v,j)  =■  I  l  M  (E.v.j;  v,j')  (3) 

v'  j'  r 

is  the  total  number  of  reactive  trajectories  without  regard  to  the  final 
vibration-rotation  state. 

The  relative  kinetic  energy  in  these  calculations  was  specified  in 
kcal/raol.  It  will  be  more  convenient  to  give  the  results  of  our  calculations 
in  terms  of  relative  velocity  rather  than  energy.  The  relation  between  velo¬ 
city  V  (in  kilometers/sec)  and  energy  E  (in  kcal/mol)  for  the  (0,H2)  system  is 


V 


2.162  E  ^ 


(4) 


The  cross  sections  computed  using  Eq.  (2)  are  plotted  versus  velocity  in  Fig. 
1.  The  solid  lines  are  least  squares  fits  of  the  computed  points  to  the  cubic 
polynomial  function 


ar(V,v,j)  -  [(a3V  +  a2)V  +  aL ]v  +  a  .  (5) 

This  fit  is  valid  in  the  range  of  velocities  VT(v,j)  <  V  <  llkm/sec,  where 
vT(v,j)  is  the  threshold  for  the  reaction.  The  coefficients  and  the  thres¬ 
holds  are  given  in  Table  2. 

B.  VIBRATIONAL  DISTRIBUTION 

The  probability  that  a  reactive  trajectory  will  terminate  in  a  particular 
vibrational  state  v"  of  the  product  OH  molecule  is  given  by 

P(v'|V,j,v)  -  l  N  (V, v, j;  v',j')/N  (V,v,j)  (6) 

r  r 
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V  Ikm/sec) 


Figure  1.  The  reactive  cross  sections  a  (V;  v,j).  The  labels  are  the  initial 
quantum  numbers  v,  j.  The  doifs  are  the  computed  values,  the  solid 
curves  are  the  least  squares  fits  using  Eq.  (5)  and  the 
coefficients  in  Table  2. 
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Table  2.  Least  squares  coefficients  and  threshold  velocities  for 

calculating  the  total  reactive  cross  sections  o  (V;  v,j)a 
by  Eq.  (5)  in  the  text. 


\) 

j 

ao 

al 

a2 

a3 

v  b 

T 

0 

0 

6.8009 

-2.5892 

0.30647 

-0.010608 

6.27 

0 

1 

8.6803 

-3.3532 

0.40406 

-0.014363 

6.27 

1 

0 

2.9638 

-1.7619 

0.30914 

-0.013505 

4.03 

1 

1 

1.6626 

-1.2653 

0.2607 

-0.011947 

3.98 

2 

The  cross  section  units  are  A  , 
Velocity  units  are  kilometers /sec. 


These  probabilities  are  plotted  versus  the  collision  velocity  in  Fig.  2.  The 
solid  lines  are  (segmented)  linear  least  squares  fits  where  the  discontinuous 
change  in  slopes  occurs  at  the  threshold  for  exciting  the  v'  -  1  state  in  Fig. 
2a  and  at  the  threshold  for  v"  -  2  in  Fig.  2b.  Above  the  thresholds  the 
probabilities  are  given  by  the  linear  fit 

P(v'|V,v,j)  -  bQ+  bx  V.  (7) 

The  coefficients  and  thresholds  are  given  in  Table  3. 

The  state-to-state  reactive  cross  sections  or(V,v,j;  v")  are  obtained  by 
multiplying  the  total  cross  sections  by  the  conditional  probability  that  if  a 
reaction  occurs  it  will  terminate  in  the  v'th  state 

°r(V,v,j;  v')  -  P(v"| V,v, j)  or(V,v,j)  (8) 

An  analytic  expression  for  this  cross  section  is  the  product  of  the  least 
squares  formulas  (5)  and  (7).  It  is 

ar(V,v,j;v')  -  (bQ+  b^)  [(a3  V  +  a2)  V  +  aj  (9) 

which  is  valid  for  velocities  in  the  range  from  the  thresholds  given  in  Table 
3  to  ~11  km/sec.  For  convenience  we  have  calculated  the  cross  sections  for 
the  case  v'  ■  I  which  are  needed  for  plume  signature  studies  and  plotted  them 
in  Fig.  3.  It  is  evident  from  this  figure  that  the  reactive  cross  section  for 
producing  OH(v'  *  1)  is  greatly  enhanced  by  a  quantum  of  vibrational  energy  in 
the  H2  reactant  molecule.  This  will  have  a  significant  effect  on  the  rate  of 
production  of  OH(v“  ■  I)  only  if  the  reactants  initially  have  a  significant 
*  I)  population,  i.e. ,  if  the  vibrational  temperature  of  the  H2  is  high 
enough. 

The  effect  of  the  initial  vibration  and  rotation  temperatures  of  the 
reactants  can  be  taken  account  of  by  defining  a  "distribution  averaged"  cross 
section. 
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OH(v')  molecule.  The  dots  are  the  computed  values,  the  solid  lines 
are  a  least  squares  fit. 


Table  3.  Lease  squares  coefficients  and 
threshold  velocities  for  calcu¬ 
lating  the  final  vibrational 
distribution  by  Eq.  (7)  in  the 
text. 


V 

j 

v' 

bo 

bl 

v_a 

T 

0 

0 

0 

1.242 

-0.0323 

7.5 

0 

0 

1 

-0.242 

0.0323 

7.5 

0 

1 

0 

1.256 

-0.0346 

7.4 

0 

1 

1 

-0.256 

0.0346 

7.4 

1 

0 

0 

0.310 

0.0003 

4.0 

1 

0 

1 

-0.827 

-0.0252 

4.0 

1 

0 

2 

-0.245 

0.0361 

6.8 

1 

1 

0 

0.117 

0.0159 

4.0 

1 

1 

1 

0.976 

-0.0431 

4.0 

1 

1 

2 

-0.237 

0.0356 

6.7 

aVelocity  units  are  kilometers/sec 


V  (km/sec) 


The  cross  section  for  the  reaction  0  +  H  (v, j)  ♦OH  (v'  ■  1)  H 
The  curves  were  calculated  using  Eq.  (8)  and  the  coefficients  in 
Tables  2  and  3.  The  labels  on  the  curves  are  v,  j. 


(10) 


o  (V, T  ,T  ;  V')  -  q“X  I  I  g  (2 j  +  1) 

v  j  J 

exp[-Ev(v)/kTv  -  Er  (y,j)/kTR]or(V,v,j;  v") 

In  this  equation  Ty  and  TR  are  the  vibration  and  rotation  temperatures,  k  is 
Boltzmann's  constant,  gj  is  the  ortho-para-statistical  weight  factor  (g^  *  1 
for  even  j,  gj  ■  3  for  odd  j),  Ey  (v)  and  ER  ( v,  j)  are  the  vibration  and 
rotation  energies  of  ^(v.j),  and  Q  is  the  normalizing  partition  function 
obtained  by  evaluating  the  sum  in  Eq.  (10)  with  the  last  cross  section  factor 
equal  to  unity.  The  vibration  and  rotation  energies  can  be  adequately 

Q 

approximated  using  the  simple  vibrating-rotor  formulas 

Ey(v)  -  mg(v  +  i  )  -  mexe(v  +  |)2  (11) 


and 


ER(v,J)  -  Bvj(j  +  1)}  Bv  *  Be  -  <x(v  +1)  (12) 

where  the  symbols  have  their  usual  meaning. 

We  have  used  Eq.  (10)  to  evaluate  the  "distribution-averaged"  cross 
section  curves  for  the  production  of  0H(v"  -  1)  for  a  fixed  rotational  tem¬ 
perature  TR  *  200°K  and  three  values  of  the  vibrational  temperature  Ty  ■  1000° 
K,  1500°  K  and  2000°  K.  These  are  plotted  versus  velocity  in  Fig.  4.  Below 
1000°K  the  contribution  of  v  *  1  to  the  reaction  rate  is  quite  small,  although 
for  velocities  less  than  the  threshold  of  7.4  km/sec  it  represents  the  only 
contribution.  Above  1000°  K  the  v  -  1  states  of  the  reactant  H2  molecules 
significantly  Increase  the  cross  sections. 

These  calculations  were  repeated  for  two  other  rotational  temperatures, 
Tr  *  100°  K  and  TR  *  300°  K,  with  very  little  change  from  the  200°  K  results 
shown  in  Fig.  4.  This  is  to  be  expected  since  the  state-to-state  cross  sec¬ 
tions  are  not  strongly  dependent  on  the  rotational  states. 

C.  ROTATIONAL  DISTRIBUTION 

In  a  given  final  vibrational  state  v",  the  probability  distribution  over 
the  rotational  states  j"  is  given  by 
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(13) 


P(riv'5  V.v.j)  -  N  (V.v.j;  v',j')  /  £  N  (V.v.j;  v'.j') 

r  y  r 

These  conditional  probabilities  were  computed  and  then  fit  to  the  analytic 
function 

P(j'|v';  V,v,j)  -  (2j"  +  1)[U  -  j'(j'  +  l)]1/2exp[-  0  j'(j'  +  l)/o]/Q 

(14) 

This  function  has  been  used  in  previous  work  to  represent  rotational 
distributions^  and  is  derived  by  information  theoretic  methods. The 
normalization  constant  Q  is  determined  by  the  condition 

{***  P( j " | v";  V, V, j)  -  1  (15) 

r-  0 

The  quantity  U  is  defined  to  be 

U  •  [E  +  Ev>j  -  Ev.  -  A]/Bv,  (16) 

where  E  is  the  relative  kinetic  energy  of  the  reactants,  E  is  the 

v,  j 

vibrational-rotational  energy  levels  of  H2,  Ey,  is  the  vibrational  energy 
levels  of  OH,  4  ■  2.79  kcal/mol  is  the  endothermicity  of  the  reaction  measured 
with  . respect  to  the  potential  energy  minima  of  the  product  and  reactant 
channels,  and  is  the  rotational  constant  of  OH  in  the  v"th  vibrational 

level.  Finally,  6  is  the  rotational  surprlsal  parameter  which  can  be  adjusted 
to  obtain  a  fit  to  the  trajectory-computed  rotational  distributions.  The 
quantity  jmx  in  Eq.  (15)  is  the  maximum  value  of  the  rotational  quantum 
number  allowed  by  energy  conservation.  It  is  the  integer  part  of  the  solution 
of  the  equation 


W-W  + »  ■ »  <17> 

The  quantity  U,  defined  by  Eq.  (16),  depends  on  the  quantum  numbers  v,j 
and  v''.  However,  the  dependence  on  j  is  weak  and  can  be  ignored.  Also,  for 
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Che  rocket  plume  signature  calculations  we  only  need  to  determine  the  rota¬ 
tional  distributions  associated  with  the  v"  -  1  vibrational  mode.  Therefore, 
we  need  only  consider  two  functional  forms  for  U,  corresponding  to  v  *  0  and 
v  -  1.  To  sufficient  accuracy  they  are 

UQ  -  (E  -  12. 5)/0. 051  (18) 


and 


UL  -  E/0.051  (19) 

where  U  and  E  are  in  units  of  kcal/mol. 

These  expressions  are  then  substituted  into  Eqs.  (14)  and  (15)  to  obtain 
the  rotational  distributions  and  To  obtain  reasonable  fits  to  the 

computed  distributions  it  was  found  that  the  parameter  9  was  dependent  only 
on  v  and  that  the  data  could  be  fit  using 

9  -  1.0,  (v  -  0)  (20) 

and 

9-0,  (v  -  1)  (21) 

Calculation  of  the  rotational  distribution  using  Eq.  (14)  is  easily  imple¬ 
mented  on  a  computer  or  even  a  programmable  calculator. 

Figure  5  shows  a  comparison  of  the  actual  trajectory-calculated  rota¬ 
tional  distribution,  plotted  as  a  histogram,  and  the  fitted  function  which  is 
shown  as  a  continuous  curve,  for  a  typical  case  v-1,  j  «  1,  v"  -  1  and  E  - 
25  kcal/mol.  Figure  6  shows  how  the  rotational  distributions  vary  as  a 
function  of  the  initial  vibration  state  v  and  the  relative  kinetic  energy  of 
the  reactants. 

The  final  partitioning  of  the  state-to-state  reactive  cross  section  into 
its  constituent  product  rotational  states  can  now  be  computed  by  multiplying 


22 


1 


Figure  5.  Product  rotational  distribution  for  the  case  v  ■  1,  j  ■  1,  v'  ■ 
and  E  •  25  k  cal/raol.  The  histogram  is  the  trajectory  results 
computed  using  Eq.  (13).  The  smooth  curve  is  the  analytic  fit 
computed  using  Eq.  (14)  with  the  surprisal  parameter  3*0.  (The 
smooth  curve  was  plotted  for  visual  convenience;  only  the  integer 
values  of  j"  have  physical  meaning.) 


Figure  6.  Produce  rotational  distributions  computed  using  the  analytic  fit 

Eqs.  (14)  -  (21).  Each  curve  is  labeled  with  the  initial  relative 
kinetic  energy  E  (keal/mol)  of  the  reactants;  v  is  the  initial 
vibrational  state. 


the  vibrationally  resolved  cross  section,  given  analytically  by  Eq.  (9),  by 
the  rotational  probability  function,  given  analytically  by  Eq.  (14),  to  obtain 

®r(V,v, j;v-, j')  -  P(j'|v';  V, v, j)  Or(V,v,j;v')  (22) 

This  formula  is  the  main  result  of  this  report.  It  represents  (implicitly) 
the  state-to-state  cross  sections  as  a  product  of  three  factors,  each  of  which 
is  a  simple  analytic  (least  squares  fit)  formula.  The  three  factors  are  the 
reaction  cross  section  given  by  Eq.  (5),  the  vibrational  probability  distribu¬ 
tion  given  by  Eq.  (7)  and  the  rotational  probability  distribution  given  by  Eq. 
(14).  This  analytic  formula,  along  with  the  parameters  given  in  Tables  2  and 
3  and  Eqs.  (18)  -  (21),  provides  an  efficient  means  for  calculating  the  state- 
to-state  reactive  cross  sections  for  the  0  +  H-,  reaction. 
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LABORATORY  OPERATIONS 


The  Laboratory  Operations  of  The  Aerospace  Corporation  Is  conducting  exper¬ 
imental  and  theoretical  investigations  necessary  for  the  evaluation  and  applica¬ 
tion  of  scientific  advances  to  new  military  6pace  systems.  Versatility  and 
flexibility  have  been  developed  to  a  high  degree  by  the  laboratory  personnel  in 
dealing  with  the  many  problems  encountered  in  the  nation’s  rapidly  developing 
space  systems.  Expertise  in  the  latest  scientific  developments  is  vital  to  the 
accomplishment  of  tasks  related  to  these  problems.  The  laboratories  that  con¬ 
tribute  to  this  research  are: 

Aerophysics  Laboratory:  Launch  vehicle  and  reentry  aerodynamics  and  heat 
transfer,  propulsion  chemistry  and  fluid  mechanics,  structural  mechanics,  flight 
dynamics;  high-tempera ture  thermomechanics ,  gas  kinetics  and  radiation;  research 
in  environmental  chemistry  and  contamination;  cw  and  pulsed  chemical  laser 
development  including  chemical  kinetics,  spectroscopy,  optical  resonators  and 
beam  pointing,  atmospheric  propagation,  laser  effects  and  countermeasures. 

Chemistry  and  Physics  Laboratory:  Atmospheric  chemical  reactions,  atmo¬ 
spheric  optics,  light  scattering,  state-specific  chemical  reactions  and  radia¬ 
tion  transport  in  rocket  plumes,  applied  laser  spectroscopy,  laser  chemistry, 
battery  electrochemistry,  space  vacuum  and  radiation  effects  on  materials,  lu¬ 
brication  and  surface  phenomena,  thermionic  emission,  photosensitive  materials 
and  detectors,  atomic  frequency  standards,  and  bioenvironmental  research  and 
monitoring. 

Electronics  Research  Laboratory:  Microelectronics,  GaAs  low-noise  and 
power  devices,  semiconductor  lasers,  electromagnetic  and  optical  propagation 
phenomena,  quantum  electronics,  laser  communications,  Ildar,  and  electro-optics; 
communication  sciences,  applied  electronics,  semiconductor  crystal  and  device 
physics,  radiometric  imaging;  millimeter-wave  and  microwave  technology. 

Information  Sciences  Research  Office:  Program  verification,  program  trans¬ 
lation,  performance-sensitive  system  design,  distributed  archi tectures  for 
spaceborne  computers,  fault-tolerant  computer  systems,  artificial  intelligence, 
and  microelectronics  applications. 

Materials  Sciences  Laboratory:  Development  of  new  materials:  metal  matrix 
composites,  polymers,  and  new  forms  of  carbon;  component  failure  analysis  and 
reliability;  fracture  mechanics  and  stress  corrosion;  evaluation  of  materials  in 
space  environment;  materials  performance  in  space  transportation  systems;  anal¬ 
ysis  of  systems  vulnerabi li ty  and  survivabi 11 ty  in  enemy-induced  environments. 

Space  Sciences  Laboratory:  Atmospheric  and  ionospheric  physics,  radiation 
from  the  atmosphere,  density  and  composition  of  the  upper  atmosphere,  aurorae 
and  airglow;  magnetospheric  physics,  cosmic  rays,  generation  and  propagation  of 
plasma*  waves  in  the  magnetosphere;  solar  physics,  infrared  astronomy;  the 
effects  of  nuclear  explosions,  magnetic  storms,  and  solar  activity  on  the 
earth's  atmosphere,  ionosphere,  and  magnetosphere;  the  effects  of  optical, 
electromagnetic,  and  particulate  radiations  in  space  on  space  systems. 


